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In the two-impurity Anderson model, the inter-impurity spin exchange interaction favors a spin 
singlet state between two impurities leading to the breakdown of the Kondo effect. We show that 
a local uniform magnetic held can delocalize the quasiparticles to restore the Kondo resonance. 
This transition is found to be continuous, accompanied by not only the divergence of the staggered 
(antiferromagnetic) susceptibility, but also the divergence of the uniform spin susceptibility. This 
may imply that the magnetic field induced quantum phase transitions in Kondo systems are in favor 
of the local critical type. 
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The study on quantum phase transitions and critical 
phenomena has been an extraordinarily active area of re- 
search in condensed matter physics and quantum field 
theory. One example which has been extensively stud- 
ied in experiments is the magnetic quantum phase tran- 
sition (QPT) in heavy fermion metals Two theo- 
retical scenarios are suggested for this QPT: the spin- 
density-wave picture based on itinerant quasiparticles, 
and the local quantum criticality mandating the local- 
ization of quasiparticles. While it is commonly acknowl- 
edged that the competition between the onsite Kondo 
coupling and the intersite spin exchange interaction, 
namely, the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction, plays the determinant role, there are few 
theoretical methods which can handle them on an equal 
footing. However, this problem is well-defined in the two- 
impurity Anderson (or Kondo) model, which presents 
such a competition effect in an exactly solvable way [H- 
[loj | . With the Kondo coupling, the impurity spin forms 
a Kondo singlet state with the spins of the conduction 
electrons and a quasiparticle resonance peak develops at 
the Fermi energy, which is described by the Kondo ef- 
fect. When the inter- impurity spin exchange interaction 
is antiferromagnetic and strong enough, the two impu- 
rity spins tend to form a singlet by themselves, against 
the formation of Kondo singlets, leading to the local- 
ization of quasiparticles. As a result, the quasiparticle 
spectra have a "pseudogap" at low energies [1, • It 
is found that, the phase transition between the Kondo 
resonance state and the inter-impurity spin singlet state 
is continuous |3|, accompanied by the divergence of the 
staggered (antiferromagnetic) spin susceptibility together 
with a discontinuous change of the spectral weight at the 
Fermi energy 0, HBl ■ 

A magnetic field has been an essential experimental 
tuning parameter to investigate the magnetic properties 
and spin correlations of condensed matter physics, which 
is also relevant to the two-impurity Anderson model. 
However, a detailed analysis on the magnetic field ef- 



fect for this model, especially close to the characteristic 
scales of the two-impurity quantum critical point (QCP), 
is still lacking, which is the purpose of this study. While 
there are some existing theoretical studies llMl3j target- 
ing the double quantum dot, they are limited to the cases 
with large magnetic fields, in which t he p hysical proper- 
ties follow the Zeeman splitting effect [ll|, [l2 1 . While it is 
known that a local staggered magnetic field can induce a 
QPT as it directly couples to the critical staggered spin 
fluctuations, it is not known the role of a local uniform 
magnetic^ field, which is usually applied in experimental 
studies Hi]. 

In this Letter, we report the first observation of a mag- 
netic field induced quantum phase transition in the two- 
impurity Anderson model from a numerical study. We 
find that a local uniform magnetic field on the two impu- 
rities can drive a transition from the inter-impurity spin 
singlet state to the Kondo resonance state, leading to the 
derealization of quasiparticles. We further show that 
this transition is continuous, accompanied by the abrupt 
change of the quasi-particle spectral weight at the Fermi 
energy and the divergence in staggered spin susceptibil- 
ity. In sharp contrast with the two-impurity QCP at zero 
field, the uniform spin susceptibility is found to be also 
divergent at this magnetic-field-induced QCP. The new 
observation is suggestive that the field-induced QCP in 
heavy fermion systems [3] does have the local nature, as 
advocated in recent QCP theories [jjj . 

The Hamiltonian for the two-impurity Anderson model 
can be written as 



H = £>4 
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which describes two interacting local orbitals fi a (Ander- 
son impurities) in hybridization with a non-interacting 
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conduction electron band Ck a with the strength Vk at 
each impurity site r^. e/ and U are the energy level 
and onsite Coulomb interaction for the local orbitals, re- 
spectively. I is a direct spin exchange interaction be- 
tween two impurities. We here consider a uniform mag- 
netic field B which acts on the impurity spins only, while 
h = g/J.BB has the dimension of energy g, [Ib are Lande 
factor and Bohr magneton, respectively. In reality, an 
applied magnetic field acts on the conduction electrons 
as well. But as along as the Zeeman energy is much 
smaller than the bandwidth, we can safely neglect this 
effect and only consider its effect on the impurity spins. 
This model has been shown 0, 0, E3] to be equivalent to 
a two-impurity two-channel model, with degrees of free- 
dom casted into the even (e) and odd (o) parity chan- 
nels. The local orbitals become / 6i0 = (fi ± /2V V2 and 
the hybridization functions for these two channels are 
r e , H - (l/2iV c )^y k 2 |e lk - r ±e- lk - r | 2 ^- ek ), where 
r = (rj. — i"2)/2. The inter-site spin exchange interaction 
can be generated by considering specific forms of Vk and 
£k (RKKY interaction) , or provided by the direct spin ex- 
change term I. We adopt the numerical renormalization 
group (NRG) [16| method with the complete-Fock-space 
NRG (CFS-NRG) method for calculations of dynam- 
ical quantities at zero temperature, including the spectral 
function Af pa = —\mGf iPa (bS), the uniform and stag- 
gered spin susceptibilities Xu, a = ((Si z ±S 2z ; Si z ±S 2z )). 
We notice that the CFS-NRG method is particularly suit- 
able for the problem with a finite magnetic field, as the 
characteristic energy scale is in the intermediate energy 
range (close to the Kondo temperature). This has also 
been evidenced from a similar study on the single impu- 
rity model [l8j . 

In our previous study on this model l(j . we have 
explicitly calculated a system with a well-defined two- 
impurity QCP (we refer the two-impurity QCP to this 
zero-field QCP in the following). Our results can be sum- 
marized as follows. We choose T et0 (u) = Tq, for which no 
RKKY interaction is generated and the single-impurity 
Kondo temperature Tr can be determined. We then add 
the direct spin exchange interaction / (to simulate RKKY 
interaction) to tune the competition between the Kondo 
effect and the inter-impurity spin exchange interaction. 
Below a critical value I c , the low energy properties are 
still due to the Kondo effect, Af(0) « l/(7iTo) and 
XUO) ~ 1/T K , but Xa(0) ~ V t *f with t *f the reduced 
(local) Fermi liquid temperature. Above I c , it is the 
inter- impurity spin singlet state with vanishing Af(Q). It 
differs from the Mott gap in the fact that there are still 
finite spectral weights at low energies: Af(ui) ~ uj 2 for 
u < Tp and a Non-Fermi liquid form for T F < ui < T s f , 
where Tp and T s f (spin fluctuation scale) correspond re- 
spectively to the two energy scales Tj, and Th identi- 
fied in our previous work [I0(. For Tq = 0.0457r£>, and 
ef = —U/2 = —D, it is found that the single-impurity 
Kondo temperature Tk = 1.0 x 10~ 3 £> and the critical 



value I c w 0.0023464D w 2.3T K . At I c , there is a sud- 
den change of the spectral weight at the Fermi energy, but 
the transition is still continuous. This is evidenced by the 
uniformly vanishing Tp ~ (I — I c ) 2 with the divergence 
of the staggered spin susceptibility xj,(0). However, the 
uniform spin susceptibility x«(0) remains finite through 
the transition. 
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FIG. 1: (color online) The scattering phase shift 5 e ^(0) as a 
function of h for different values of /. The phase shifts at the 
Fermi energy for different channels and spins are the same 
due to the symmetries A £rr (u}) = A oct (cj) and A e ^(uu) = 
A e i(—cj). The data labelled with "Kondo" are subtracted 
from Costi's paper on the single-impurity Kondo model [2(j. 
The inset shows the critical value h c as a function of / — I c for 



/ > Ic. Two dotted lines are h c = A.5Tk(I/Ic 
and he = 6.0T K {I/Ic - 1) 2/3 (black). 



,1/2 



(red) 



We add a local uniform magnetic field h to the above 
system to examine its effects. In Fig. [1] we show the 
results of the scattering phase shift 5 pa (Q) at the Fermi 
energy, determined from A/(0) = (1/ttTo) sin 2 S pa (Q), as 
functions of h for various values of /. For 1 = 0, this is 
equivalent to a single-impurity Kondo problem and in- 
deed our results are in agreements with those obtained 
for the single-impurity Kondo model [20(. The scatter- 
ing phase shift by the exact Bethe-Ansantz method is 
5 h (0) ~ n /2-h for h<.T K while ^(0) - l/\og{h/T K ) 
for h ^> Tk- The latter relation is in agreement for all fi- 
nite Is but with 1/ \og(h/T s f). For / < I c , Sh(0) is always 
finite and is enhanced from the single-impurity case for 
the same h. For / > I c , 6h(0) vanishes when h — 0, char- 
acterizing the inter-impurity spin singlet state. A finite 
but small h induces a small Sh (0) or a small quasiparticle 
weight at the Fermi energy: it is found that 5h(0) ~ h 
as verified by a log-log plot (not shown). When h is in- 
creased to a critical value h c , we observe a sudden jump 
of the phase shift from a tiny value to a large value of 
the order unity, which indicates a transition rather than 
a crossover between the inter-impurity spin singlet state 
and the Kondo resonance state. The relation between h c 
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and I — I c is shown in the inset of Fig. Q] We find that 
h c ~ (I - I c ) 1/2 for (I — I c )/I c <C 1 and a noticeable 
deviation for (I — I c )/I c > 0.1. Such a deviation is also 
identified in T£ ~ (I - I c ) a when h = [h|. 



h/T K =1.2 
h/T K =1.2814 
h/T K =1 .4 
h/T K =6 




l/(7rro) but smaller. This is similar to the two-impurity 
QCP, except here the jump amplitude is smaller. The 
staggered spin susceptibility has the same behavior as 
well: it becomes divergent when h — > h c . The significant 
difference lies in the uniform spin susceptibility: it is also 
divergent at h c in this field-induced QCP compared with 
finite Xu(Q) ~ in the two- impurity QCP. While we 

can subtract the energy scales from l/[4x„ >a (0)] as deter- 
mining Tj(, a more reliable method to determine the low 
energy scale is from a scaling analysis, which is shown 
in Fig. for x"( w )- O nce the energy is scaled with a 
certain scale Th for different h, the low energy part of 
falls into a universal curve. Similar behavior can 
be observed in but its high energy part does not 

appear to scale (or not universal). We further plot the 
obtained Th as a function of \h — h c \: Th can be fitted as 
T h ~ \h — h c \ 2 for both h < h c and h > h c . Th — > as 
well as the divergence in both the staggered and uniform 
spin susceptibilities are clear evidences that this magnetic 
field induced QPT is continuous. 
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FIG. 2: (color online) Spectral functions Af(uu) (a) and the 
imaginary parts of the uniform and staggered spin suscepti- 
bilities (b) as functions of energy for different values of h in 
the inter- impurity singlet regime, I — 3.07V. In (a), the solid 
and dotted lines represent the spin up and spin down degrees 
of freedom, respectively (the even and odd channels are the 
same). In (b), they respectively represent the uniform and 
staggered spin susceptibilities. The inset shows the real parts 
of the spin susceptibilities at zero energy (static) as functions 
of h. The inverse values are shown for convenience to be 
compared with the Kondo scale. 

To show that it is a continuous phase transition, we 
present a detailed analysis for the I = 3.0Tk case. In 
Fig. [2j we show the results of the spectral functions 
Af tPa (u]), and the uniform and staggered spin suscep- 
tibilities Xu,a{&)- As the parity symmetry is not broken 
with a uniform magnetic field h, Af_ ea {ui) — Af^ oa (uj). 
For a small magnetic field, for instance h — 1.2Tk, the 
spectral weight at the Fermi energy is only slightly en- 
hanced. Around T s f, where spin fluctuations reach max- 
imum, the peak positions (on the positive energy range) 
of the spectral functions for the spin-up and spin-down 
degrees of freedom are different, with the energy differ- 
ence roughly given by h. For a large magnetic field, for 
instance h — QTk, which is bigger than either the Kondo 
scale or the RKKY scale, this is similar to the single- 
impurity Kondo model in the presence of a magnetic 
field. Spin-up and spin-down resonance peaks are lo- 
cated respectively at positive and negative energies, and 
the energy difference (or the gap) is 2h, which is the hall- 
mark of Kondo resonance (for a non-interacting orbital, 
the gap from Zeeman splitting is h). Around a critical 
value h c ss 1.2814TR- , the spectral weight at the Fermi en- 
ergy jumps from a tiny value to a value comparable with 




FIG. 3: (color online) Scaling behavior of the staggered spin 
susceptibility (dotted lines) for various values of h near h c . 
The rescaled uniform spin susceptibility (solid lines) is also 
shown. From the scaling, Th can be obtained and is plotted 
in the inset as a function of \h — h c \/h c - The line is a fitting 
Th/Tk = 0.28|/i/7ic - 1| 2 . 



When I I c , h c vanishes and the field- induced QCP 
merges with the two-impurity QCP at zero-field. We 
then need to understand why the divergence in \ u van- 
ishes. The results for the spectral functions and the spin 
susceptibilities for / = 2.3Tk are shown in Fig. 01 In- 
deed, any small h induces the full Kondo resonance at 
the Fermi energy with the finite Fermi temperature Th, 
which is fitted as Th ~ h 4 . Wc also observe enhancement 
of uniform spin fluctuations, which is manifested as a flat 
part above Th, x"( w ) ~ Ch- @h increases as h increases. 
When Th vanishes as h — > 0, Ch also vanishes. As a re- 
sult, Xu is not divergent. However, for any / > I c , as Ch 
remains finite, \u indeed diverges as Th vanishes when 
h — > h c . The variation of Ch also explains that X«( w ) 
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does not scale for oj > Th [cf. Fig. 3]. 
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FIG. 4: (color online) Spectral functions Af(uu) (a) and the 
imaginary parts of the uniform and staggered spin susceptibil- 
ities (b) as functions of energy for different values of h in the 
two- impurity quantum critical regime, I = 2.3Tk- The line 
representations in (a) and (b) are the same as in Fig. 2. The 
inset shows Th as a function of h, which is obtained the same 
way as in Fig. 3. The line is a fitting T h /T K = 0.18{h/T K ) 4 . 

The continuous QPT induced by a local uniform mag- 
netic field in the two-impurity Anderson model was not 
predicted by either the conformal field theory Q or the 



bosonization construction 0, 
channel QCP where T h ~ h 2 



9H- Compared with the two- 
2l| , the analogy of the mag- 



netic field in the two-channel Kondo impurity model is 
the staggered magnetic field in the two-impurity model, 
i.e., h s (Si z — S2z), which directly couples to the critical 
staggered spin fluctuations. Our results suggest that the 
uniform magnetic field h also couples effectively to the 
staggered spin fluctuations, as evidenced from the diver- 
gence of Xa- We can make the following statements on 
this coupling term based on our results. 1) The control 
parameter is modified as / — I c — ah 2 , to be consistent 
with the exponents deduced from our numerical data in 
different regimes. In other words, the field-induced QCP 
is the same in nature as the two-impurity QCP at zero 
field. From another perspective, for a given finite h, tun- 
ing / can also lead to a QPT but the critical value of I c 
is shifted up. 2) It involves the uniform spin fluctuations 
as h is directly coupled to. This also accounts for the 
divergence of the uniform spin susceptibility. In other 
words, this divergence is induced rather than the driving 
mechanism. Similar divergences are also observed in the 
non-local superconducting fluctuations and the current 
fluctuations between two impurities. 3) It vanishes or 
becomes irrelevant as h vanishes, as necessary to explain 
the loss of divergence in Xu at the two-impurity QCP. 

In a lattice generalization from the self-consistently 
solved two-site cluster, the uniform and staggered spin 



susceptibilities in this cluster correspond to the lattice 
spin susceptibility at momentum points Q — (0, 0, 0) 
and Q n = (tt, tt, tt) for a three-dimensional (3D) lattice. 
The divergence in Xa relates to the antiferromagnetic in- 
stability near Q^. However, in 3D, the spin density of 
states near vanishes, unlike in 2D. The local spin 
susceptibility, which is a sum of contributions from all 
momentum points, is not divergent. This may relate to 
the spin-density-wave type of transitions. If the transi- 
tion is driven by the magnetic field, we learn from this 
study that the divergence at Q„ can induce the diver- 
gence in Qq as well. The local spin susceptibility in this 
case diverges even for 3D. This corresponds to the local 
critical type of transitions. This might be related to the 
local quantum behaviors observed in YbRh 2 Si2 tuned by 
the magnetic field 22 1. 
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